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In recent years, tremendous progress has been 
made in exploring and exploiting the analogy of 
classical light and matter waves for fundamen- 
tal investigations and applications [Ij. Extend- 
ing this analogy to quantum matter wave optics 
is promoted by the intrinsic non-linearity of in- 
teracting particles, a stepping stone towards non- 
classical states [Ills]. In light optics, twin-photon 
beams [4j are a key element in providing the non- 
local correlations and entanglement required for 
applications like precision metrology and quan- 
tum communication [5j. Similar sources for mas- 
sive particles have so far been limited by the 
multi-mode character of the involved processes 
or a predominant background signal [6-13J. Here 
we present highly efficient emission of twin-atom 
beams into a single transversal mode of a waveg- 
uide potential. The source is a one-dimensional 
degenerate Bose gas [14j in the first radially ex- 
cited state. We directly measure a suppression 
of fiuctuations in the atom number diff*erence be- 
tween the beams to 0.37(3) with respect to the 
classical expectation, equivalent to 0.11(2) after 
correcting for detection noise. Our results under- 
line the high potential of ultracold atomic gases as 
sources for quantum matter wave optics and will 
enable the implementation of schemes previously 
unattainable with massive particles |El I15j - [T9J . 

Binary cohisions between atoms provide a natural 
means to generate dual number states of intrinsically cor- 
related atoms [T5^. Experimental schemes include spon- 
taneous emission of atom pairs by collisional deexcita- 
tion or four- wave mixing [10]-[12]. Stimulated emis- 
sion into twin-modes has been demonstrated in seeded 
four- wave mixing [H |6] , and parametric amplification in 
optical lattices [71 [8] or spinor condensates J3l IT9H2T] . 
For motional states, suppression of relative number fluc- 
tuations could so far only be demonstrated for multi- 
mode twin-atoms [12]. A different route to non-classical 
states is provided by ensembles in multi-well potentials, 
that become number-squeezed during their time evolu- 
tion [221 Eg. 

Here, we demonstrate how collisional deexcitation of a 
one-dimensional degenerate Bose gas can be used to effi- 
ciently create matter wave beams of twin- atoms. The re- 
stricted geometry of a waveguide potential enforces emis- 




FIG. 1. Schematic of the excitation and emission process, (a) 
The quasi-BEC is transferred from the ground state |0,0;0) 
into |1,0;0), the first excited state of the trapping potential 
along the radial ^/-direction. This is accomplished by means of 
fast non-adiabatic movement of the potential minimum along 
an optimized trajectory (inset). The excited state decays by 
emission of twin-atoms into the radial ground state modes 
|0,0;±/co). (b) After excitation and pair emission, the cloud 
is released from the trapping potential and imaged after ex- 
pansion. The central part of the system clearly shows the 
spatial structure of the radially excited state (blue). Two 
clouds containing the twin- atoms (red) are emitted. 



sion of the beams into a single transversal mode, in anal- 
ogy to an optical parametric amplifier [4 . We prepare 
the initial population inversion to a radially excited state 
by shaking the trap following an optimal control strategy. 
Time-of-flight fluorescence imaging is used to directly ob- 
serve the suppressed relative number fluctuations in the 
emitted beams. 

The starting point of our investigations is a dilute, 
quantum degenerate gas of neutral ^^Rubidium atoms 
magnetically trapped in a tight waveguide potential with 
a shallow axial harmonic confinement {vx = 16.3 Hz) 
on an atom chip [24]. Our scheme relies on an effec- 
tive two-level system in the radial vibrational eigenstates 
of the waveguide. This is accomplished by creating un- 
equal level spacings in the radial z-plane by radio fre- 
quency dressing [25j, which introduces anharmonicity 
and anisotropy. The resulting single-particle first and 
second excited state energies are E^^l = /i- [1.83, 2.58] kHz 
and Ey'^J = h • [3.82, 5.22] kHz. Due to the increasing 
level spacings, the ground state \ny^nz]kx) = |0,0;0) 
{riy^z and kx denoting the radial quantum numbers and 
the axial momentum, respectively) and the first excited 
state along |1,0;0) have the lowest energy difference 
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FIG. 2. Dynamics of the excitation and emission process: 
comparison between theory and experiment, (a) Optimized 
control trajectory for the trap movement along y. (b) Mea- 
sured momentum distribution. Fluorescence images inte- 
grated along kx over a region encompassing the central cloud. 
Average over 5 experimental runs. (c) Fit of the time- 
dependent momentum distribution, scaled to the experimen- 
tal total density at each time step, (d) Calculated momentum 
distribution, using a one-dimensional GPE model. As the 
model does not take into account the pair emission process, 
agreement to the experiment is expected approximately up to 
the end of the excitation pulse, (e) Red line: Population of 
the emitted clouds obtained from the same data set as (b). 
Red crosses: population of the emitted clouds obtained from 
a separate measurement with 100 experimental runs each (er- 
ror bars are the ensemble standard deviation). Black line: 
theoretical estimation for spontaneous processes only. 



among all possible combinations (figure [T^), establishing 
a closed two-level system. 

Using standard techniques we generate a Bose gas of 
typically 700 atoms at a temperature T < 40 nK ^ 
h/ks- 830 Hz (obtained independently from fits to a non- 
excited degenerate gas [26 and its residual thermal frac- 
tion [27]). The thermal occupation of state |1,0;0) is 
negligible, and the chemical potential, quantifying the 
mean- field interaction is /i ^ h ■ 500 Hz <C Ey^"^ . Our Kb 
sample is therefore a one-dimensional, weakly interacting 
quasi-Bose-Einstein condensate (quasi-BEC, interaction 
parameter 7 ~ 0.008 [14 ), which can be seen as a lo- 
cally coherent matter wave [28 with a coherence length 
approximately one order of magnitude below the system 
size. 

Having prepared the gas, we create a population in- 
version by transferring the quasi-BEC almost entirely to 
state |1,0;0) (see figure [l^). The transition is driven by 
shaking the trap along the radial ^/-direction on the scale 



of the ground state size (~ 100 nm). The trajectory (total 
duration 5 ms, see figure [2|i) has been optimized employ- 
ing an iterative optimal control algorithm (see methods) . 
In the experiment, the displacement is achieved by driv- 
ing a current in an auxiliary chip wire, parallel to the 
main trapping wire. 

We monitor the radial momentum distribution of the 
quasi-BEC by releasing the cloud from the trapping po- 
tential at different times t during and after the exci- 
tation pulse. Images are taken after 46 ms of ballistic 
expansion (figure [2]3), using a single-atom-sensitive flu- 
orescence imaging system [29]. After the excitation, we 
observe a small residual beating between the macroscopi- 
cally occupied |1, 0; 0) state and a remaining non-excited 
population in |0,0;0). From a flt to the beating pat- 
tern (flgure ^) we estimate an efficiency of the coher- 
ent transfer of rje ~ 97% and deduce the energy differ- 
ence e = h • 1.78kHz between |0,0;0) and |1,0;0) (see 
methods). A more detailed, quantitative model of the 
excitation dynamics is given in [30 . The slight devia- 
tion between Ey^^ and e is explained by particle interac- 
tions. A calculation based on the one-dimensional Gross- 
Pitaevskii-Equation (GPE) (figure [2]i) shows excellent 
agreement to the observed dynamics. 

The population inversion to |1, 0; 0) represents a highly 
non-equilibrium state of the system, analogous to a laser 
gain medium after a pump pulse. For the ensuing re- 
laxation, the only allowed channel is a two-particle col- 
lisional process, emitting atom pairs with opposite mo- 
menta. In contrast to experiments with free-space col- 
lisions [12] or two-dimensional gases [9 , the con- 
stricted geometry and non-degenerate level scheme of 
our source restricts the outgoing matter waves to the ra- 
dial ground state of the waveguide, yielding twin-atom 
beams in a single transversal mode. Within a binary 
collision, two atoms are scattered from |1, 0; 0) |1, 0; 0) 
to |0, 0; +/co) IO7 0; — ^0)7 where energy conservation re- 
quires the final momenta to be centred around ±/co = 
±V2rne/h. The emission process can be understood as a 
matter wave analogue to a degenerate optical paramet- 
ric amplifier, where the initially empty twin- modes are 
seeded by vacuum fluctuations [4] and gain an exponen- 
tially growing population if phase matching conditions 
are fulfllled. Due to flnite size, the axial multi-mode char- 
acter of the quasi-BEC source [14], and the depletion of 
the initial state, a quantitative description of the process 
is challenging. In ref. [30^ , a model based on a density ma- 
trix expansion has been developped, which mostly covers 
those effects. A comparison of the observed rates to a 
simple calculation using Fermi's golden rule (see meth- 
ods) is shown in flgure |2^ and demonstrates the insuffl- 
ciency of a purely spontaneous model (in contrast to the 
flndings of [9 for a transversal multi-mode system). 

Once the trap potential is switched off and the atoms 
propagate freely (flgure [l]3), the twin-beam modes can 
be detected essentially background-free in fluorescence 
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FIG. 3. Atom cloud image analysis, (a) Typical experimental 
image of ^ 700 atoms released from the trap 7 ms after start- 
ing the excitation sequence. The cloud is allowed to expand 
for 46 ms, making the initial momentum distribution acces- 
sible. The quasi-BEC in the excited state |1,0;0) is clearly 
distinct from the emitted clouds at momenta zthko. Units 
are photons per pixel. The blue box indicates the integra- 
tion range for the data shown in figure (b) Average over 
^ 1500 images similar to (a). The colour scale is logarith- 
mic (dB referenced to peak density). The regions used for 
correlation analysis are indicated as red boxes, (c) Normal- 
ized, radial momentum distributions of the central (blue) and 
emitted (red) clouds. Average of 50 images of clouds released 
at t = 6 ms. As comparison, the distribution of a non-excited 
cloud is shown (black, average over 100 images), (d) Normal- 
ized profile of (b) along kx (red dots) and three-peak fit (black 
line) based on stochastic simulations ^26^ . 
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FIG. 4. Correlation analysis, (a) Histogram of observed signal 
imbalances s between the emitted clouds, in units of the bino- 
mial standard deviation atin = {pSy^'^. The curves indicate 
normal distributions corresponding to the experimental result 
of = 0.11(2) (black, solid), the limits of perfect correlation, 
where only detection noise remains (red, solid), of uncorre- 
lated signals, defining the reference point for (blue, solid), 
and a binomial distribution for pS trials (black, dashed), (b) 
Observed signal imbalance variances for data bins correspond- 
ing to different total signal in the emitted clouds S. Error 
bars are the standard error. The lines correspond to those in 
panel (a). The corrected variances are given by the vertical 
distances between the data points and the detection noise. 



images (figure |3^,b) as they separate from the source. 
In figure [Sj), the radial momentum distribution of the 
twin-beams is compared to an independent measurement 
of the initial |0,0;0) cloud. The small deviation is at- 
tributed to a slight overlap of the central cloud into the 
integration regions of the emitted clouds (red boxes in 
figure [sJd) and interactions with the the mean field of 
the quasi-BEC in the |1,0;0) state. Furthermore, an ex- 
cited cloud at t = 6 ms is shown. From the axial position 
of the side peaks (figure [sji), we can deduce an emitted 
atom momentum of /cq = 27r • 0.883(3) pm~^, equivalent 
to e = h • 1.78(1) kHz, in perfect agreement with the 
value determined from the beating fit (figure ^) . The 
width of the emitted clouds (figure [3]i) is increased by 
a factor of ~ 1.4 with respect to the source cloud. At 
momenta corresponding to e' ^ h ■ 3.9 kHz, very weakly 
populated additional atom clouds (< 1 atom per image) 
are observed on an averaged picture (figurejsfD). They im- 
ply a transient population of |2,0;0) during the control 
sequence, which directly decays into the radial ground 
state. In the following analysis, they are merged with 
the atoms originating from |1,0;0). 



The non-classical correlation in the emitted twin-atom 
beams is revealed by a sub-binomial distribution of the 
number imbalance n = Ni — N2 between atoms detected 
at ±/co. The variance of n can be expressed as = ^^TV, 
where N denotes the mean total atom number in the 
emitted clouds. The noise reduction factor quantifies 
the suppression of with respect to a binomial distri- 
bution, and thus the amount of correlation between the 
populations A^i and N2. 

In the fluorescence images, we count photons in re- 
gions encompassing the emitted clouds, which have been 
released at t = 7 ms. For given atom numbers A^i,2, the 
expectation values for the photon numbers are 5'i^2 = 
pNi^2 + where p = 12.3(9) denotes the average num- 
ber of photons per atom and h/2 accounts for background 
events. Our main observable is the variance cri of the sig- 
nal imbalance s = Si — 82- Its expectation value for a 
binomial distribution of atoms is given by a^^^ = pS^ 
where S = §1 ^ S2 — h. From the experimental data as 
shown in figure |4^ we obtain an uncorrected reduction 
factor crl/al.^ = 0.37(3). However, a significant contri- 
bution to does not originate from the atom number 
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fluctuations, but from the detection process itself. For 
fluorescence imaging as employed in our experiment, this 
contribution can directly be calculated from photon shot 
noise and detection background (see methods). It is ac- 
counted for by subtracting a correction from . From 
the corrected variances, we infer a reduction factor of 
= (cr^ — cy'^)/crlin = 0.11(2), which is the main result 
of this paper. It is equivalent to an intensity squeezing 
in the sense of [4]. In strong contrast to the suppressed 
relative fluctuations, applying an analogous calculation 
on the variance <j| of the summed signal in the emit- 
ted clouds S = Si -\- S2 (binned into groups with similar 
total atom number) yields super-Poissonian fluctuations 
(cr| — a'^)/pS ~ 7, again highlighting the presence of 
bosonic ampliflcation. 

To study the correlation data in more detail, we bin the 
experimental shots according to the emitted atom signal 
S and calculate the variances and mean signals S for 
each bin (consisting of typically 100 runs) separately, as 
shown in flgure ^p. The differences between the data 
points and the detection noise represent the corrected 
variances as introduced above. They appear to be inde- 
pendent of aS, which is supported by x^-test results on 
various plausible models. As any uncorrelated emissions 
should scale with this suggests that the non-zero value 
of can be explained by a slight additional background 
signal, e.g. due to the residual overlap of the excited 
quasi-BEC and the emitted clouds. 

The availability of single-mode twin-atom beams adds 
an essential building block for quantum matter wave op- 
tics. As our scheme does not rely on the internal struc- 
ture of the atoms, it can be applied to any sufficiently 
controllable system of interacting bosons. Possible appli- 
cations include interferometry with dual-Fock states [15 , 
Hong-Ou-Mandel type experiments [17] or continuous- 
variable entanglement [5]|19]. Inclusion of internal (e.g. 
hyperfine states with appropriate scattering properties) 
or few-mode external (e.g. two-mode double well states) 
degrees of freedom seems a viable strategy to generate 
non-local entangled states of massive particles for Bell- 
type measurements [16} [TSl [T9l [311 [32]. 
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FIG. 5. Schematic of the atom chip layout. The waveguide 
potential is formed by the current Itrap through the main 
trapping wire (black) and a static magnetic field By. On a 
separate chip layer, currents lax in broad wires (green) pro- 
vide axial confinement. An external field Bx completes the 
loffe-Pritchard configuration. The radio frequency dressing 
currents Irf are applied to wires (red) in parallel to the trap- 
ping wire. Finally, the modulation of the trap position is 
accomplished by a current Imod in an auxiliary wire (blue). 
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METHODS 

Preparation of the trapping potential 

Both the optimized excitation scheme needed for transfer- 
ring the quasi-BEC into |1,0;0) (see next section), and re- 
stricting the emission to a single transverse mode |0,0;±/co) 
require a sufficiently anharmonic trapping potential along y 
with increasing level spacings. Hence, the initially radially 
symmetric loffe-Pritchard field configuration created by our 
chip wire configuration (see [33] and figure [5| is being mod- 
ified by radio frequency dressing [251 . Typically used for 
creating double well potentials, this technique also allows for 
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the introduction of anharmonicity and anisotropy to a sin- 
gle trap when the dressing strength is kept slightly below the 
point where actual splitting of the potential occurs. We ap- 
ply an ac current of Irf = 23 mA peak-to-peak amplitude at 
a detuning of ^ = —54 kHz with respect to the atomic Larmor 
frequency near the trap minimum (uq — 824 kHz) to two wires 
running on each side of the main trapping wire at a distance 
of 55 pm. 

The resulting potential can be calculated numerically by 
means of a Floquet analysis [35] . In the two radial directions it 
can be approximated by quartic polynomials of the form E = 
P4r^ -\-p2r^- In the ^-direction, along which the excitation is 
performed the coefficients are p4 = /i- 13.1 Hz/ro^ and p2 — h- 
343Hz/ro^. In the 2;-direction perpendicular to the excitation 
motion the coefficients are pa — h ■ 10.4Hz/ro^ and p2 — 
h ■ 793Hz/ro^ Here, ro = 172 nm is the mean radial ground 
state radius as calculated by the Gross-Pitaevskii equation. 

Along the axial x-axis, the trap frequency given in the text 
is determined by observation of a deliberately excited sloshing 
mode of the quasi-BEC. 



Optimized excitation of the condensate 

To transfer the cloud into the radially excited state |1, 0; 0) 
we displace the radial trap minimum along an optimized tra- 
jectory. 

This movement is achieved by applying a current of typi- 
cally less than 10 mA to a wire parallel to the trapping wire 
at 140 pm distance (see figure[5|. The resulting magnetic field 
is mostly oriented along the z-direction at the cloud position, 
moving the potential predominantly along y "24^. The cal- 
culated offset is 26nm/mA along y and 9nm/mA along 
where the contribution along z does not significantly distort 
the excitation process (see below). 

The most efficient way of moving the trap for excited state 
preparation is obtained from numerical calculations which 
employ optimal control 36 of the Gross-Pitaevskii equation. 
Optimal control theory is a powerful tool which allows to 
minimize a given cost functional with the constraint that the 
system is governed by the corresponding equations of motion. 
In our case, we iteratively solved an optimal control system 
determined from a Lagrangian framework [361 137j . For de- 
sired state trapping one maximizes the overlap of the wave 
function at the final time with a given desired state. This has 
been used before in theoretical works about transferring or 
splitting a BEC by continuously transforming a trap poten- 
tial from an initial to a final shape, without exciting the BEC 
[38 , 39 . In the present work, we employ the same techniques, 
but choose the first excited state of the Gross-Pitaevskii equa- 
tion at the final trap position as desired state. 

For such excited state preparation, it is crucial to use an 
anharmonic trap, as for a weakly interacting system in a har- 
monic trap, displacements generate displaced ground states, 
i.e., coherent states. The present optimal control calculations 
have been performed in ID (along y) for simplicity, using an 
effective ID interaction parameter [40]. This works well, since 
the dynamics in the axial x-direction is orders of magnitudes 
slower than in the radial ones. Moreover, the potential is 
sufficiently anisotropic in the radial plane, such that the z 
direction is not significantly affected by the movement, even 
though it is not strictly performed along y (see above). We 
performed 2D simulations of the Gross-Pitaevskii equation 



and found only small deviations as compared to the ID case. 

To estimate the excitation efficiency from experimental 
data we fit a model of the time-dependent momentum den- 
sity along y in the main cloud (i.e. not including the emitted 
atoms) after the excitation pulse (see figure 2b, c in the main 
text), according to 

n{ky,t) = (1 - r]d{t))\y^l - r]eipo{ky) + ^ry^e*^/^'Vi(^?/)|^ 
+ r]4t)\Mky)\\ 

where the total density of the experimental data is nor- 
malized for each time step. The first line corresponds to 
the coherent two-level dynamics between the radially excited 
state ^l^i{ky) — (/cy|l,0;0) and the ground state ipo{ky) — 
(/cy |0, 0; 0) at a beating frequency given by the single-particle 
energy difference e ~ The efficiency of excitation is 

denoted as rje. In the second line, an incoherent admixture 
of the ground state is added with a time-dependent popula- 
tion rfdit)^ which corresponds to dephased atoms having un- 
dergone emission into the |0,0;±A:o) modes and subsequent 
collisions with the excited main cloud. The momentum- space 
wave functions V^i and V^o correspond to eigenstates of the 
one-dimensional Gross-Pitaevskii equation at a typical atom 
number. From the fit we deduce an efficiency of the coherent 
excitation of ?7e = 0.97 and a dephased fraction linearly rising 
from ?7cz(5.2ms) = 0.08 to ?7cz(8ms) = 0.28. For the energy 
difference we obtain e — h- 1.78 kHz. 



Emission dynamics model 

We compare the fraction of emitted high-momentum atoms 
to a simple model, taking into account spontaneous emis- 
sion only, and assuming an atom number corresponding to 
the average of the shown dataset [Ntot ~ 700). To ob- 
tain the theory curve, we solve rate equations between the 
states |0,0;0), |1,0;0) and |0,0;it/co). The excitation rate 
from |0,0;0) to |1,0;0) cannot be directly obtained from 
GPE calculations, as both the populations of the radial states 
as well as their wave functions are time-dependent. It is 
thus approximated by a constant rate during the excitation 
pulse (see ref. 30 ). On the other hand, to treat the emis- 
sion from |1,0;0) to |0,0;±A:o), we assume an axial den- 
sity distribution according to ^Ij and employ Fermi's golden 
rule for indistinguishable bosons to obtain a weakly atom- 
number-dependent two-body decay constant of the order of 
rem[A^eccc] 0.05s~^. This allows to estimate the emission 
rate as iV = 2Tem[Nexc(t)]Nexc(tY ^ where Nexc denotes the 
number of atoms in |1,0;0), calculated as described before. 
While during the initial phase, where Nexc,N <^ Ntot the 
model works well, it fails at later times where Nexc, N < Ntot, 
indicating that stimulated processes have to be taken into 
account, analogous to parametric amplification in quantum 
optics. 



Fluorescence detection and noise corrections 

To detect the atoms after release from the trapping po- 
tentials we employ a fluorescence detector introduced in 2^ . 
The expanding cloud falls through a thin horizontal light sheet 
(vertical waist radius wq = 20 pm) after 46 ms of expansion 
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time. Light scattered by the atoms is collected by an objective 
outside the vacuum vessel and imaged on an EMCCD camera. 
The extraordinarily low background of the system enables 
single-atom sensitivity while maintaining the dynamic range 
needed for imaging of dense Bose-Einstein condensates. For 
the excitation power employed (which is of the order of the 
atomic saturation intensity at the centre of the light sheet), 
p = 12.3di0.9 photons per atom are detected on average. This 
number is deduced from cross-calibration by comparison of a 
sufficient number of calibrated absorption imaging and fluo- 
rescence shots of clouds prepared under identical conditions. 
It is compatible with independent estimates obtained from 
physical properties of clouds released from a nearly isotropic 
trap [27 . Slight dependencies of p on the spatial position 
within the image due to CCD etaloning and inhomogeneity 
of the illumination can be corrected by using a reference im- 
age. 

For exact analysis of the pair correlation as described be- 
low, also the variance of the distribution of photons per atom 
cTp has to be estimated. The first contribution to this vari- 
ance is the photon shot noise cr|^ = p. Furthermore, the 
excess noise due to the stochastic amplification in the EM- 
CCD detector has to be taken into account and increases the 
detection noise by a factor of two ^amp — 2p, in ac- 

cordance to independent characterization measurements us- 
ing Poissonian light sources and the value for the detection 
background (see below). Finally, a further broadening of the 
distribution of p is caused by the diffusion of atoms within the 
light sheet [29]. However, this contribution is hard to obtain 
from experimental results (as we cannot prepare single atoms 
deterministically as e.g. in [J^) and can only be estimated 
from simulations. To avoid overestimation of detection noise, 
which would spuriously reduce , we assume cr^ = cr^mp and 
set (Tp/p — 2 in our analysis. 

Apart from the distribution of p, another detection contri- 
bution to originates from the residual background signal 
b. We extract b and al from regions directly adjacent to the 
main analysis regions, which contain the emitted clouds and 
scale the obtained values accordingly. We obtain al/b — 2.14, 
which can be understood as combined effect of shot and am- 
plification noise and residual readout noise. Again, to avoid 
over-correction, we make sure that the background treatment 
does not lead to spurious reduction of when the analysis 
regions are enlarged artificially. 



Correlation analysis 

As a first step to determine the amount of population cor- 
relation between the emitted clouds, we count the photon 
numbers 5'i,*S'2 (indicating atom numbers A^i, A^2) within re- 
gions encompassing them (see figure 3b in the main article). 
Choosing the integration range along the axial x-direction is a 
sensitive task, as a too small region will not capture all emit- 
ted atoms, deteriorating the result. On the other hand, one 
has to take into consideration residual overlap of the emitted 
clouds with the quasi-BEC. As expected, we find that setting 
the limit to the position of the density minimum between 
main and emitted clouds minimizes the obtained value of 
Still the residual overlap appears to be one of the limiting 
factors to the measurable correlation. The outer limits are 
chosen in a way that the weak clouds emitted from |2,0;0) 
(on average less then one atom per image) are still contained 



within the analysis regions. Those overlap significantly with 
the main emitted clouds, thus limiting the ability to assess 
the correlations between them separately. 

To calculate from the accessible quantities al and as 
defined in the text, we have to estimate the total contribution 
of detection noise cr^ from the quantities cr^ and (Tp as defined 
in the previous section. We can decompose into a popu- 
lation fluctuation term and a detection term using the law of 
total variance: 

Var(s) =Var(E(s|n)) + E(Var(s|n)) 

=Var(E(p)n) + E(Var(p)Ar + Var(6)), 

where b and p are the background signal and the number of 
photons per atom as discussed above, s = 5*1 — & , n = A^i — 
N2, N — N1+N2. E(X) and Var(X) denote expectation value 
and variance of the random variable X. For our analysis, 
assummg Var(X) = cr\ and E(X) = X, this reads 

(jl^P^fN + alN + al 
^pfS + vS + al 

2 I 2 
— ^ C^bin^ C^d^ 

and thus cr^ = vS + cr^. Solving for leads to the formula 
used for the analysis. 

To estimate the uncertainty of we propagate the stan- 
dard error of the shot-to-shot mean of b^a^^S and p, the lat- 
ter obtained from an independent measurement as described 
earlier. The error of a1 — {k — Sm=i('^^ ~ ^)^' where 
k denotes the number of experimental runs, is calculated as 
a?.(2/(fc-l))i/^ 
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